Mon. Not. R. Astron. Soc. 000, [|-| (2000) Printed 1 February 2008 (MN WT$L style file vl.4) 



Global VLBI Observations of Compact Radio Sources in 
M82 

A. R. McDonald, 1 * T. W. B. Muxlow, 1 A. Pedlar, 1 M. A. Garrett, 2 
K. A. Wills, 3 S. T. Garrington, 1 P. J. Diamond, 1 P. N. Wilkinson 1 

1 University of Manchester, Jodrell Bank Observatory, Macclesfield, Cheshire SK11 9DL. 

2 Joint Institute for VLBI in Europe, Postbus 2, 7990 Dwingeloo, The Netherlands 

3 Department of Physics, University of Sheffield, Sheffield S3 7RH. 



Accepted ?. Received ? 



ABSTRACT 

Observations of the starburst galaxy, M82, have been made with a 20-station global 
VLBI array at A18cm. Maps are presented of the brightest young supernova remnants 
(SNR) in M82 and the wide-field mapping techniques used in making images over a 
field of view of ~1 arcminute with 3 milliarcsecond resolution are discussed. A limit 
has been placed on the power law deceleration of the young SNR, 43.31+592 with an 
index greater than 0.73 ±0.11 from observations with the European VLBI Network. 
Using the global array we have resolved compact knots of radio emission in the source 
which, with future global observations, will enable better constraints to be placed on 
the expansion parameters of this SNR. 

The latest global observations have also provided high resolution images of the 
most compact radio source in M82, 41.95+575. We determine an upper limit to the 
radial expansion rate along the major axis of 2000 km s _1 . However, the new images 
also show structure resembling that of collimated ejection which brings into question 
the previous explanation of the source as being due to the confinement of a supernova 
by a high density circumstellar medium. 

It is apparent that we are now able to image the brightest supernova remnants 
in M82 with a linear scale which allows direct comparison with galactic SNR such as 
Cassiopeia A. 

Key words: galaxies : starburst - galaxies : individual : M82 - radio continuum : 
galaxies - techniques: interferometric - techniques : image processing 



1 INTRODUCTION 

M82 is generally considered to be the archetypal 'starburst' 
galaxy. Far-infrared (FIR) and radio luminosities im ply a 



high star-formation rate and at a distance of 3.2 Mpc (Bur- 
bidge, Burbudge & Rubin 1964) it lends itself well to study 



with radio interferometry on a parsec scale. The radio stud- 
ies of M82 have revealed a population of compact radio 
sources (Unger et al., 1984; Kronberg et al., 1985) which 
were tentatively identified as being either radio supernovae 
(RSn) or supernova remnants (SNR). Subsequent observa- 
tions with MERLIN at a resolution of 50 mas (0.75 pc) re- 
solved most of the compact sources and the parsec-scale 
shell-like structure observed in many of the object s con- 
firmed them to be young SNR (Muxlow et al. 1994). Fur- 



thermore, recent EVN observations at 15 mas (0.2 pc) res- 



* email: amd@jb.man.ac.uk 



olution (Pedlar et al. 199£) have further resolved the most 
compact sources and by comparison with a 1986 epoch an 
expansion speed of ~9500 kms -1 has been measured for the 
compact shell-like remnant, 43.31+592. Obviously it would 
be preferable to measure the size evolution of this remnant 
on much shorter timescales than the ~10 years required by 
the EVN observations. A global VLBI array at A18cm pro- 
vides a maximum resolution of ~3 mas and enables these 
measurements to be performed on timescales of ~2 years 
and also provides an opportunity to determine the deceler- 
ation of the remnants as they interact with the interstellar 
medium (ISM). The SNR in M82 and other nearby galaxies 
provide a means of understanding the star-formation pro- 
cesses in a starburst environment and also enable the ISM 
in a starburst to be probed on parsec scales. 

Since the SNR are produced by massive stars with life- 
times of order 10 6 years, then these SNR trace the recent 
star formation history of a galaxy. If the supernova rate can 
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Table 1. Positions used for the three correlator runs on the global 
VLBI data. 





Source 


Rig 


ht Ascension 


Declination 








(J2000) 


(J2000) 


1 


41.95+575 


09^ 


55 m 50.694 s 


69° 40' 43.69" 


2 


43.31+592 


Q9 h 


55 m 52.030 s 


69° 40' 45.42" 


3 


44.01+596 


09 h 


55 m 52.720 s 


69° 40' 45.78" 



be established, then by assuming an initial mass function it 
is possible to calculate a total star formation rate to com- 
pare with that derived from other methods, such as the to- 
tal far-infrared or non-thermal radio luminosity. Previous 
VLBI observations of M82 have concentrated on the bright- 
est and most compact source, 41.95+575 (Bartel et al., 1987; 
Trotman 1996). However, we are now able to process the 
large amounts of data needed to image a large field of view. 
Therefore, it is now possible to image a field of view of ~1 
arcminute with the 3 mas resolution which the global array 
provides at A18cm. 



2 THE OBSERVATIONS AND IMAGE 
PROCESSING 

M82 and a nearby phase calibrator, 0955+697, were ob- 
served at 1652 MHz with 20 antennas of a global VLBI 
network on 28th November, 1998. The network comprised 
all 10 antennas of the Very Long Baseline Array (VLBA), 
7 antennas of the European VLBI Network (EVN), the two 
NASA Deep Space Network (DSN) antennas at Robledo, 
Spain and Goldstone, USA and a single VLA dish. The tele- 
scopes received both left and right hand circular polarisa- 
tions (LCP & RCP) which, after correlation, were combined 
to form the total intensity maps presented in this paper. 
The observations were made with the VLBA/MklV record- 
ing system and a total bandwidth of 16 MHz. The observa- 
tions were carried out in phase reference mode with a cycle 
of 8 minutes on the target (M82) and 3 minutes on the phase 
calibrator. At the time of correlation, this was the largest 
VLBI experiment ever performed and the 190 baselines pro- 
vided by the combinations of 20 stations resulted in excellent 
UV-coverage which in turn provided a greater image fidelity 
than had previously been achieved. 

The data were correlated at the National Radio Astron- 
omy Observatory, Socorro, New Mexico, using the VLBA 
processor. The correlator used a pre-averaging time of 2 
seconds and a frequency resolution of 0.25 MHz. Three sep- 
arate correlations of the data were performed in order to 
help provide the ~1 arcminute field of view required to im- 
age the SNR. Each correlation was centred on the position 
of a bright compact source and these positions are given in 
Table 1. 

All subsequent processing was carried out using the 
NRAO AIPS package. The visibility amplitudes were cal- 
ibrated using the system temperature and gain information 
provided for each telescope and the residual delays, fringe 
rates and antenna gains were determined from 0955+697 
and interpolated and applied to M82. In order to avoid 
time-averaging and bandwidth smearing, the data were kept 



in the form of 64 continuous 0.25 MHz channels and were 
not averaged in time. Initially, the data correlated at the 
position of the brightest source, 41.95+575, were Fourier 
transformed and CLEANed using the AIPS task IMAGR. 
The solar activity at the time of observation was high 
so that the phase stability during the observing run was 
not good. Therefore, the phase solutions derived from the 
nearby phase calibrator, 0955+697, were not sufficient to 
completely remove atmospheric phase fluctuations in the di- 
rection of M82. Therefore, further improvements were made 
using standard self-calibration techniques. The delay, fringe 
rate and phase solutions derived from this initial fringe- 
fitting and self-calibration procedure were then applied to 
the other two pointings which were then Fourier transformed 
and CLEANed. The increase in resolution from 15 mas with 



the EVN (Pedlar et al. 1999) to 3 mas results in a decrease 
in surface brightness sensitivity. Hence, only two of the five 
sources detected with the EVN are detected at full resolu- 
tion. However, by applying a Gaussian taper to the data 
in the UV-plane, maps were made at 10 mas and 15 mas 
resolution for comparison with the EVN images. 

The 3 mas resolution maps made from the data cor- 
related at the position of the brightest source, 41.95+575, 
had a rms noise off-source of 45 /iJy beam -1 . However, the 
same phase and amplitude corrections when applied to the 
other two correlations were only sufficient to produce maps 
with noise levels of 68 piJy beam -1 at the same resolution. 
This difference in sensitivity between different correlations 
illustrates the importance of having a bright source with 
which self-calibration can be performed near to the pointing 
centre of the observations. Therefore, maps made from the 
data correlated at the positions of 43.31+592 and 44.01+596 
were only used as a consistency check. However, the ability 
to make wide-field maps and hence image 43.31+592 from 
data correlated at the position of 41.95+575 emphasises the 
power of correlating the data with a wide-field mode (i.e. 
short integration times and small channel bandwidths). A 
comparison between the 3 mas resolution maps of 41.95+575 
made from pointings 1 and 2 is shown in figure [j]. To check 
for consistency at slightly lower resolution, figure |^ shows 
a comparison between the maps of 43.31+592 also made 
from pointings 1 and 2 and restored with a CLEAN beam 
of FWHM 4 mas. 

By applying a purely uniform weighting sc heme (with 
the Brigg's robustness parameter (Briggs 1995) set to -2) 
to the UV data, a maximum resolution of ~2.5 mas was 
achieved with a corresponding decrease in sensitivity to 
88 fiJy beam . However, of all the sources in M82, only 
41.95+575 has a sufficient brightness to be detectable at this 
resolution and sensitivity and the resulting map is shown in 
figure |^( a). In addition, figure ^(b) shows the same source 
at slightly lower resolution. 

A weighting scheme intermediate between the uniform 
weighting required to image 41.95+575 and pure natural 
weighting was applied to the data and images with a sensi- 
tivity sufficient to image the shell-like remnant, 43.31+592 
were made. Figure ^ shows the image of 43.31+592 at a 
resolution of 4 mas. 



3 DISCUSSION 
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Figure 1. 18cm global VLBI images of the compact source 41.95+575 in M82 at 3 mas resolution, (a) was made from data correlated 
to the position of 41.95+575, (b) was made from data correlated to the position of 43.31+592 and (c) is the difference between (a) and 
(b). The greyscale ranges from 0.2 to 2 mjy beam -1 and the contour levels are (-1, 1, 2,..., 9, 10) X 0.2 mjy beam -1 . 
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Figure 2. 18cm global VLBI images of the young supernova remnant, 43.31+592 in M82 at 4 mas resolution, (a) was made from data 
correlated to the position of 41.95+575, (b) was made from data correlated to the position of 43.31+592 and (c) is the difference between 
(a) and (b). The greyscale ranges from 0.05 to 0.3 mjy beam -1 and the contour levels are (-1, 1, 1.5, 2, 2.5, 3) X 0.1 mjy beam . 



3.1 The nature of the compact source, 41.95+575 

3.1.1 Radio morphology 

The brightest, most compact source in M82, 41.95+575 has 
been extensively studied wi th VLBI for around fifteen years. 
Observations at 2.3 GHz (|Bartel et al. 198/j) and 5 GHz 



(Wilkinson & de Bruyn 199C) showed structure which was 
interpreted to be an elongated shell and accounting f or the 
fact that 41.95+575 was detected at epoch 1965.6 (Bash 
196S ) , limits were placed on its expansion speed of < 14000 



km s . An extensive study of the source has been under- 
taken at A18cm by Trotman (1996) and by comparing global 
VLBI images at epochs 1987.9 and 1990.7 a tentative expan- 
sion speed between the two 'hotspots' of 2500 ± 550 km s -1 
was measured. 

The new images of 41.95+575 show significant structure 
on smaller scales than have previously been observed with 
VLBI. In addition to the well established separation between 



the north-eastern and south-western 'hotspots', significant 
structure is now observed linking the two hotspots as shown 
in figure ^. The most compact structure shown in figure ^ 
was suggested in the 5 GHz EVN image of Wilkinson & de 
Bruyn (1990) and has now been confirmed. From the 5 GHz 
map of Wilkinson & de Bruyn (1990) we estimate a sep- 
aration between the brightest components of 19 ± 1 mas. 
This compares with a similar estimate from the latest maps 
of 22 ± 1 mas. Therefore in the 13.1 years between these 
epochs, the source has expanded by 0.05 ± 0.02 pc (using 
a distance to M82 of 3.2 Mpc). However, given the uncer- 
tainties in identifying equivalent source structure between 
epochs and frequency bands, we prefer to express this as an 
upper limit to the expansion rate of 4000 km s -1 (a radial 
expansion rate of 2000 km s -1 ). 

Various models have been proposed to account for the 
unusual morphology of 41.95+575. Blondin et al. (1996) sug- 
gested that an axisymmetric density distribution in the wind 
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Figure 3. Images of the most compact radio source in M82, 41.95+575. (a) was made from uniformly weighted data with Brigg's 
robustness set to -2 and restored with a 2.7 X 2.3 mas elliptical beam at a position angle of -45°. The greyscale ranges from 0.2 to 1 
mjy beam -1 and the contour levels are (-1, 1, 1.5, 2, 2.5, 3, 3.5, 4) X 0.3 mjy beam -1 , (b) was made using an intermediate weighting 
scheme (robustness = 0) and restored with a 3.1 X 2.4 mas elliptical beam at a position angle of -44°. The greyscale ranges from 0.1 to 
1.5 mjy beam -1 and the contour levels are (-1, 1, 2,..., 9, 10) X 0.15 mjy beam -1 . 
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Figure 4. Image of 43.31+592 made with uniformly weighted data (robustness = 0) and restored with a circular 4x4 mas 2 beam. The 
greyscale ranges from 0.1 to 0.3 mjy beam -1 and the contour levels are (-1, 1, 1.5, 2, 2.5, 3) X 0.1 mjy beam -1 . Labels A to D indicate 
the four most significant knots of radio emission and the cross in the centre marks the average position of these knots. The diamond 
shapes represent the expected positions of these knots 10 years after the 1998.9 epoch assuming a constant self-similar radial expansion 
speed of 10000 km s -1 . 
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from the progenitor star can lead to a protrusion emerging 
along the symmetry axis. They cite 41.95+575 as an exam- 
ple of this based on the VLBI images of Bartel (1997) and 
Wilkinson & de Bruyn (1990). However, other models which 
predict nonspherical morphologies include that of Khokhlov 
et al. (1999) which involves the triggering of a supernova 
by supersonic jets at the centre of the progenitor star. The 
high resolution of the new global VLBI image reveals struc- 
ture which may be interpreted as collimated outflow and we 
speculate that this may represent a jet-like structure such 
as that generated by the model of Khokhlov et al. (1999). 
However, it is not yet clear how this model could explain 
the high degree of confinement which the source appears to 
be experiencing. 



3.1.2 Flux density variability 

Of all the compact sources in M82, 41.95+575 has shown the 
highest degree of variability over the last ~25 years. The to- 
tal flux density of the source has been consistently monitored 
with both the VLA and MERLIN and it has been found that 
at both 5 GHz (C-band) and 1.6 GHz (L-band) the flux 
density has bee n decaying at a constant rate of ~8.5% yr _1 
(Protman 1996). Based on the total flux density as mea- 
sured by these instruments, the predicted flux density for 
41.95+575 at epoch 1998.9 is ~55 mjy, taking into account 
the 8.5% yr _1 decay rate. However, the total flux density 
recovered from the numerous cycles of self-calibration was 
43.8 ± 2.3 mjy which only represents ~80% of the predicted 
flux density. However, we do not attribute this to a larger 
than expected drop in its flux density; it is more likely due 
to missing extended flux on the long VLBI baselines. More 
recent MERLIN observations (epoch 1999.1) at 5 GHz have 
shown no evidence for any departures from the 8.5% yr _1 
decay rate (result not presented here). 



3.1.3 Summary of results on 41-95+575 

The simplest explanation for the appearance of 41.95+575 
remains that of a supernova event taking place within a high 
density molecular cloud which provides the high level of con- 
finement and results in the low expansion speeds which we 
derive. The high ambient density of the medium surround- 
ing the source may also account for the X-ray luminosity, 
since Stevens et al. (1999) propose the the radio source is 
coincident with a point source as observed by ROSAT. The 
higher resolution of the Chandra X-Ray Observatory will 
make this identification more robust once an accurate regis- 
tration with a radio image has been performed. 

It is clear that an adequate supernova-related model 
which fully explains both the radio morphology and flux- 
density variability of this source has not yet been con- 
structed. Therefore, we must also ask the question as to 
whether the source is really supernova-related. 



3.2 The evolution of the shell-like remnant, 
43.31+592 

3.2.1 Limits on deceleration from low-resolution maps 

The global VLBI map of 43.31+592 represents the third 
epoch for which a map has been produced of sufficient res- 



olution to image the structure of the source. In an earlier 
paper, Pedlar et al. (1999) described the results of the first 
two epochs which consisted of observations conducted by the 
European VLBI Network in 1986 and 1997. By taking slices 
at different position angles across the shell they measured 
angular diameters of 28.5 ± 1.0 mas and 42.1 ± 1.5 mas for 
epochs 1986.9 and 1997.5 respectively. These angular sizes 
correspond to linear sizes of 0.442 ± 0.016 pc and 0.653 ± 
0.023 pc at 3.2 Mpc. Using the nomenclature of Huang et 
al. (1994) we parameterise the size evolution of the remnant 
as, 



kT 



(1) 



where D pc is the shell diameter in parsecs, T yr is the age 
of the remnant in years, k is a constant and 8 is a 'deceler- 
ation parameter'. Since the date of the supernova event is 
unknown, there are three free parameters to which we need 
to fit our data. Therefore, the current data are not sufficient 
to constrain the parameters. In addition, the uncertainty in 
the age of the remnant does not allow a correction for the 
positive bias in the size measurement of the shell due to 
the finite VLBI beamsize as discussed by Marcaide et al. 
(1997). Figure | shows four example size evolution curves 
which would be consistent with our data. However, since 
the remnant clearly existed in 1972, as is shown in the 8.1 
GHz map of Kronberg & Wilkinson (1975), it can be seen 
that the smaller values of 8 are not allowed by our measure- 
ments. 

A more exhaustive analysis using only two measured 
diameters may be conducted by deriving the relationship 
between the age of the remnant at the 1986.9 epoch and the 
deceleration parameter as 



AT 



„ m ln{D 2 /Dl) 



(2) 



where AT is the epoch separation in years and D2/D1 is 
the ratio of the remnant sizes at the two epochs. Therefore, 
the minimum value of 8 assuming a supernova event taking 
place before 1972 can be calculated to be 0.73 ± 0.11. 

This result emphasises the importance of making im- 
ages with higher spatial and hence linear resolution. Even 
with observations spaced by ~10 years (around a third of 
the remnant's inferred age) it is not possible to more accu- 
rately constrain the remnant's evolutionary parameters with 
linear resolutions like those produced by the EVN or VLBA 
alone. 



3.2.2 Prospects for proper motion studies from global 
VLBI maps 

Figure § shows the 1998.9 epoch global VLBI map of 
43.31+592 at a resolution of 4 mas. Labels A to D indicate 
the four most significant 'knots' of radio emission which we 
detect. Two-dimensional Gaussian fits were performed to 
these knots and the results are indicated by the size and 
position angle of the crosses. The parameters derived by 
the Gaussian fits to the components are shown in table ^| 
However, we do not claim that these knots are intrinsically 
Gaussian in nature and this method has only been used so 
that an indication of the expected future evolution of the 
remnant can be provided. The diamond shapes in figure ^ 
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Figure 5. Size evolution plots for the compact young supernova remnant, 43.31+592. The curves show four possibilities for the evolution 
with deceleration parameters from 0.4 to 1.0 (free expansion). 



Table 2. The results of 2D Gaussian fits to the brightest compact knots of radio emission in the young supernova remnant, 43.31+592. 
The right ascensions and declinations are quoted relative to 09' 1 55 m 00 s +69°40'00". It should be noted that the absolute astrometric 
precision of these positions has been degraded since self-calibration has been used. However, the relative precision has been retained. 



Label 


Right Ascension 


Declination 


Major Axis 


Minor Axis 


Position Angle 




(J2000) 


(J2000) 


(mas) 


(mas) 


(degrees) 


A 


52.03183 ± 0.00007 


45.4157 ± 0.0007 


13.3 ± 1.7 


6.5 ± 0.8 


15.5 ± 6.6 


B 


52.02988 ± 0.00009 


45.4099 ± 0.0008 


11.0 ± 1.8 


7.0 ± 1.1 


168 ± 14 


C 


52.02692 ± 0.00019 


45.4506 ± 0.0008 


15.1 ± 2.6 


8.8 ± 1.5 


124 ± 12 


D 


52.02542 ± 0.00010 


45.4431 ± 0.0012 


17.6 ± 2.8 


7.2 ± 1.2 


11.4 ± 6.4 



indicate the positions of these knots ten years after the ini- 
tial global observations, assuming a self-sim ilar expansion 
speed of 10000 km s _1 ( |Pedlar et al. 199S| ) . However, two 
years should be sufficient to detect the proper motions of 
these components since 10000 km s _1 corresponds to 0.66 
mas yr _1 at a distance of 3.2 Mpc. A more exhaustive anal- 
ysis of the high resolution images will be performed once a 
second epoch map at the same resolution has been obtained. 



3.2.3 Comparison with other sources 

Figure ^ shows a V LA image of the galactic su pernova rem- 
nant, Cassiopeia A ( Braun, Gull fc Perley 1987 ) on the same 



linear scale as the compact sources in M82. Also shown 
schematically is the young radio supernova, SN 1993J in 
M81. The size of the circle corresponds to the size of SN 



1993 J as measured by Marcaide et al. (1997) in October 
1996. This figure emphasises the difference between studying 
galactic supernova remnants and those in M82. Cassiopeia 
A has an age of around 320 years ( Ashworth 198C ) compared 
with an age of ~30 years for 43.31+592. It is therefore clear 
that studies of the young supernova remnants in M82 can 
fill an important gap in the studies of the evolution of su- 
pernova remnants between the very young (e.g. SN1993J) 
and those with ages of a few hundred years. However, any 
differences that are due to the M82 supernovae evolving in 
a starburst ISM may also become apparent. 
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43.31+592 



Figure 6. Comparison of a VLA image of Cassiopeia A (Braun, Gull & Perley 1987), the youngest supernova remnant in the Milky Way, 
with 43.31+592, the youngest supernova remnant in M82. Both SNR are shown on the same linear scale. Also shown on the same scale 
is the compact radio source, 41.95+575 and the young radio supernova SN 1993J in M81. The scale has been calculated by assuming a 
distance to Cas A of 2.8 kpc and a distance of 3.2 Mpc to M82 and M81. The greyscale for Cas A ranges from to 0.01 Jy pixel -1 , the 
greyscale for 43.31+592 ranges from 0.07 to 0.2 mjy beam -1 and the greyscale for 41.95+575 ranges from 0.1 mjy beam -1 to 1.5 mjy 
beam -1 . 



3.3 Other sources 

Of the remaining three sources detected by the EVN at 
epoch 1997.5 only 44.01+596 and 45.17+612 were reliably 
detected in this expermiment. 



be observed with any confidence. As with the other sources, 
more observations are required at later epochs to determine 
any expansion of the 'shell' and hence aid its identification 
as a supernova remnant or weak AGN. 



3.3.1 44.01+596 

Of all the compact radio sources in M82, 44.01+596 is the 
most lik ely candidate as the host of a weak active galactic 
nucleus (Wills et al. 1999). However, when imaged with the 
EVN at the 1997.5 epoch it showed structure more consis- 
tent with an identification as a supernova remnant. Wills 
et al. (1999) argued that this does not exclude the AGN 
possibility since Sgr A* also shows similar radio continuum 
structure at a similar linear scale. Figure ^ shows a compar- 
ison between the 1997.5 epoch EVN image of 44.01+596 (15 
mas resolution) with the 1998.9 epoch global VLBI image 
(10 mas resolution). The same structure that is seen in the 
lower resolution image is clearly present at higher resolu- 
tion, although the signal to noise ratio in the 10 mas map is 
not sufficient that any apparent changes in morphology can 



4 CONCLUSIONS 

These global VLBI observations have produced the highest 
resolution images yet made of the young supernova remnants 
formed in the nuclear starburst within M82. At the full res- 
olution of ~3 mas images have been made of the brightest 
and what are inferred to be the youngest sources. For the 
compact young supernova remnant, 43.31+592, we draw the 
following conclusions. 

(1) The measurements of the shell diameter from 15 



mas resolution EVN images (Pedlar et al. 1999) places a 



limit on the deceleration parameter. Since the remnant is 
present in the 1972 image of Kronberg & Wilkinson (1975), 
a limit of 8 > 0.73 ± 0.11 has been calculated. As yet, this 
limit does not allow for a physical analysis of the supernova 
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Figure 7. Comparison between the 1997.5 epoch EVN image of 44.01+596 (15 mas resolution) with the 1998.9 epoch global VLBI image 
(10 mas resolution). The contour levels for both images are (-1, 1, 2,..., 9, 10) X 0.1 mjy beam . 



model, but with future observations tighter constraints will 
be placed on the deceleration of this source. 

(2) The higher resolution images presented in this paper 
will enable proper motion studies to be carried out on the 
compact knots of radio emission as shown in figure [|. Further 
global VLBI observations are scheduled to take place in late 
2000, two years after the initial observations. This represents 
the minimum time required to detect significant proper mo- 
tion of the knots assuming a radial expansion speed of 10000 
km s _1 . 

(3) For the first time we have been able to image the 
supernova remnants in M82 on a similar linear scale to that 
which we can image the galactic supernova remnant, Cas- 
siopeia A. 

In addition, we draw the following conclusions for the 
compact radio source, 41.95+575. 

(1) The new image of 41.95+575 reveals structures 
which represent deviations from its previous description as 
an elongated shell and we speculate that these may either 
represent collimated outflow as predicted by the jet-induced 
explosion models of Khokhlov et al. (1999) or a prominence 
of the type predicted by Blondin et al. (1996). 

(2) We calculate an upper limit to the radial expansion 
rate along the major axis of 2000 km s" 1 . 

The young supernova remnants in M82 fill an important 
gap in the evolutionary timescale of supernova remnants be- 
tween the very young (e.g. SN1993J) and those with ages 
measured in hundreds of years (e.g. those in the Milky Way 
and Magellanic clouds). Therefore, M82 will continue to act 
as an important laboratory for the study of young SNR and 
any differences in evolution which are due to the starburst 
ISM may become apparent. 
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